Next Article in Journal
Monotone Positive Solutions for Nonlinear Fractional Differential Equations with a Disturbance Parameter on the Infinite Interval
Previous Article in Journal
Usefulness of the Ordinal Logistic Biplot: Analysis of the Path Taken towards a Circular Primary Sector in Spain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

(R, S)-(Skew) Symmetric Solutions to Matrix Equation AXB = C over Quaternions

1
School of Computer Science and Engineering, Faculty of Innovation Engineering, Macau University of Science and Technology, Avenida Wai Long, TaiPa, Macau 999078, China
2
Macau Institute of Systems Engineering, Faculty of Innovation Engineering, Macau University of Science and Technology, Avenida Wai Long, TaiPa, Macau 999078, China
3
School of Mathematics and Data Science, Minjiang University, Fujian 350108, China
4
Department of Mathematics, University of Manitoba, Winnipeg, MB R3T 2N2, Canada
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(2), 323; https://doi.org/10.3390/math12020323
Submission received: 29 November 2023 / Revised: 11 January 2024 / Accepted: 16 January 2024 / Published: 18 January 2024
(This article belongs to the Special Issue New Results in Matrix Analysis and Applications)

Abstract

:
( R , S )-(skew) symmetric matrices have numerous applications in civil engineering, information theory, numerical analysis, etc. In this paper, we deal with the ( R , S )-(skew) symmetric solutions to the quaternion matrix equation A X B = C . We use a real representation A τ to obtain the necessary and sufficient conditions for A X B = C to have ( R , S )-(skew) symmetric solutions and derive the solutions when it is consistent. We also derive the least-squares ( R , S )-(skew) symmetric solution to the above matrix equation.

1. Introduction

A real quaternion (also called Hamilton quaternion) is usually expressed as
H = a 0 + a 1 i + a 2 j + a 3 k | i 2 = j 2 = k 2 = i j k = 1 , a 0 , a 1 , a 2 , a 3 R .
It is a four-dimensional division algebra over the real number field R . During the past decades, matrices over quaternions have played important roles in signal processing, aerospace, color image processing, and many other areas [1,2,3,4,5,6].
Recall that an involutory matrix R is a square matrix satisfying R 2 = I , an identity matrix. Clearly, ± I are trivial involutory matrices. A matrix A C m × n is said to be an ( R , S ) -symmetric (resp. ( R , S ) -skew symmetric) matrix if there exist nontrivial involution matrices R and S such that R A S = A (resp. R A S = A ) (see [7,8]).
( R , S ) -(skew) symmetric matrices are widely used in linear system theory, numerical analysis, and physics [9,10,11,12,13,14]. Trench [8] discussed the ( R , S ) -(skew) symmetric solutions to the complex matrix equation X B = C . Dehghan and Hajarian [7] also considered the ( R , S ) -(skew) symmetric solutions to some complex matrix equations. They derived the solvability conditions for the existence of solutions and obtained the general solutions when the matrix equations are solvable. Zhang and Wang [15] derived the ( R , S )-(skew) symmetric extreme rank solutions to the quaternion system A X = B , X C = D . Some special cases of ( R , S )-(skew) symmetric matrices, such as (anti-)centrosymmetric matrices, P-(skew) symmetric matrices, generalized reflexive matrix, and reflexive (antireflexive) matrices, have been discussed [13,14,16,17,18]. For instance, Zhou et al. [18] studied the problems of centrosymmetric matrices over C , where R , S are the J, which are on the secondary diagonal and zeros elsewhere. Wang et al. [17] considered the P-(skew) symmetric solutions to a pair of quaternion matrix equations.
The matrix equation A X B = C is of substantial research significance due to its wide range of applications and its relevance in solving fundamental problems in various disciplines. The ( R , S )-(skew) symmetric matrices include many kinds of important special matrices, like centrosymmetric matrices, P-(skew) symmetric, generalized reflexive matrix, and reflexive (antireflexive) matrices; these kinds of solutions have far-reaching implications in areas ranging from mathematics and engineering to computer science and data analysis [13,14,16,17,18].
Motivated by the work and research significance mentioned above, we consider the ( R , S )-(skew) symmetric solutions and least-squares ( R , S )-(skew) symmetric solutions to the quaternion matrix equation
A X B = C ,
where A H p × m , B H n × q , C H p × q are known matrices and X H m × n is an unknown matrix.
The paper is organized as follows: In Section 2, we first introduce some preliminary results. Then we derive the solvability conditions of the equation A X B = C for ( R , S ) -(skew) symmetric solutions. In Section 3, the least-squares ( R , S ) -(skew) symmetric solutions are given. Finally, we provide a numerical example in Section 4.
Throughout this paper, we propose some notations. For a matrix A, A T , A , and r ( A ) denote the transpose, conjugate transpose, and rank of A separately. Moreover, A stands for the Moore–Penrose inverse of A. I n will be the n × n identity matrix.

2. The Solvability

The real representation method is one of the standard and efficient ways to solve questions over quaternions. There are several real representations; see, for example, [19,20,21]. In this paper, we will use the following X τ :
For X = X 1 + X 2 i + X 3 j + X 4 k H m × n , X i R m × n , i = 1 , 2 , 3 , 4 ,
X τ = X 1 X 2 X 3 X 4 X 2 X 1 X 4 X 3 X 3 X 4 X 1 X 2 X 4 X 3 X 2 X 1 .
It is easy to verify that this real representation can convert an ( R , S )-(skew) symmetric matrix into a real ( R τ , S τ )-(skew) symmetric matrix. For further discussions of our problem, we introduce the following orthogonal matrices:
Q n = 0 I n 0 0 I n 0 0 0 0 0 0 I n 0 0 I n 0 , G n = 0 0 I n 0 0 0 0 I n I n 0 0 0 0 I n 0 0 ,
T n = 0 0 0 I n 0 0 I n 0 0 I n 0 0 I n 0 0 0 .
Now, we summarize some properties of the above real representation in the following lemma:
Lemma 1.
Let A , B H m × n , C H n × s , a R . Then
(a) 
( A + B ) τ = A τ + B τ , ( a A ) τ = a A τ ;
(b) 
( A C ) τ = A τ C τ ;
(c) 
Q m T A τ Q n = A τ , G m T A τ G n = A τ , T m T A τ T n = A τ ;
(d) 
A τ = A τ T , A 1 τ = A τ 1 ;
(e) 
Let A H m × m . Then A τ commutes with Q m , G m , T m .
Definition 1.
Assume that R H m × m and S H n × n are nontrivial involutory matrices. A H m × n is an ( R , S ) -symmetric (resp. ( R , S ) -skew symmetric) matrix when A satisfies R A S = A (resp. R A S = A ).
Now, we are in a position to derive some necessary and sufficient conditions for the matrix Equation (1) to have ( R , S )-symmetric (resp. ( R , S ) -skew symmetric) solutions and provide the solutions when the matrix equation is consistent. Let R H m × m and S H n × n be nontrivial involutory matrices. By (d) of Lemma 1,
R τ = ( R 1 ) τ = ( R τ ) 1 ± I ,
S τ = ( S 1 ) τ = ( S τ ) 1 ± I ;
thus, R τ , S τ are also nontrivial involutory matrices. For R τ and S τ , according to [8], we can find positive numbers r , k and matrices P C 4 m × r , Q C 4 m × ( 4 m r ) , U C 4 n × k , V C 4 n × ( 4 n k ) such that
P P = I , Q Q = I , R τ P = P , R τ Q = Q ,
and
U U = I , V V = I , S τ U = U , S τ V = V .
Next, we denote
U ^ = U ( I + S τ ) 2 , V ^ = V ( I S τ ) 2 , P ^ = P ( I + R τ ) 2 , Q ^ = Q ( I R τ ) 2 .
For A τ R 4 p × 4 m and B τ R 4 n × 4 q of (1), we perform the following decomposition:
A τ P Q = A 1 A 2 , A 1 C 4 p × r , A 2 C 4 p × ( 4 m r ) ,
U ^ V ^ B τ = B 1 T B 2 T T , B 1 C k × 4 q , B 2 C ( 4 n k ) × 4 q .
Now, we have our main results as follows:
Theorem 1.
Let A H p × m , B H n × q , C H p × q . Then, there are three equivalent statements:
(a) 
The matrix Equation (1) has an ( R , S )-symmetric solution X H m × n ;
(b) 
The matrix equation
A τ Y B τ = C τ
has an ( R τ , S τ )-symmetric solution Y R 4 m × 4 n ;
(c) 
The following rank equalities hold:
r A 1 A 2 C τ = r A 1 A 2 , r B 1 B 2 C τ = r B 1 B 2 ,
r C τ A 1 B 2 0 = r 0 A 1 B 2 0 , r C τ A 2 B 1 0 = r 0 A 2 B 1 0 .
Furthermore, if the matrix Equation (1) is solvable, then
X = 1 16 I m i I m j I m k I m T Y + Q m Y Q n T + G m Y G n T + T m Y T n T I n i I n j I n k I n
is an ( R , S )-symmetric solution to (1), where
Y = P Q Y P U 0 0 Y Q V U ^ V ^
with
Y P U = Y ˜ P U + S 1 F J U E H T 1 + F A 1 V 1 + V 2 E B 1 ,
Y Q V = Y ˜ Q V + S 2 F J U E H T 2 + F A 2 W 1 + W 2 E B 2 ,
S 1 = I r 0 , S 2 = 0 I 4 m r , T 1 = I k 0 ,
T 2 = 0 I 4 n k , J = A 1 A 2 , H = B 1 B 2 ;
U , V 1 , V 2 , W 1 , W 2 are arbitrary, and Y ˜ P U , Y ˜ Q V are particular solutions to the matrix equation  A 1 Y P U B 1 + A 2 Y Q V B 2 = C τ .
Proof. 
Let R H m × m and S H n × n . For the matrix Equation (1), we convert it into the matrix equation A τ Y B τ = C τ by the real representation method. First, it can be shown that each ( R τ , S τ )-symmetric solution Y of the real matrix Equation (3) can generate an ( R , S )-symmetric solution X of the original matrix Equation (1).
Suppose (3) has an ( R τ , S τ )-symmetric solution Y, i.e., R τ Y S τ = Y . Y may not have the structure of a real representation; thus, we need to construct Y with the structure of a real representation from Y. According to Lemma 1, we have the following three equations:
Q p T A τ Q m Y Q n T B τ Q q = Q p T C τ Q q ,
G p T A τ G m Y G n T B τ G q = G p T C τ G q ,
T p T A τ T m Y T n T B τ T q = T p T C τ T q .
Since Q p , Q n , G p , G n , T p , T n are nonsingular,
A τ ( Q m Y Q n T ) B τ = C τ ,
A τ ( G m Y G n T ) B τ = C τ ,
A τ ( T m Y T n T ) B τ = C τ ,
and by Lemma 1(e), we have
R τ Q m Y Q n T S τ = Q m R τ Y S τ Q n T = Q m Y Q n T ,
R τ G m Y G n T S τ = G m R τ Y S τ G n T = G m Y G n T ,
R τ T m Y T n T S τ = T m R τ Y S τ T n T = T m Y T n T .
Therefore, Q m Y Q n T , G m Y G n T , and T m Y T n T are also the ( R τ , S τ )-symmetric solutions of (3), and so is
Y = 1 4 Y + Q m Y Q n T + G m Y G n T + T m Y T n T .
Suppose that Y can be written in the form of a block matrix:
Y = Y 11 Y 12 Y 13 Y 14 Y 21 Y 22 Y 23 Y 24 Y 31 Y 32 Y 33 Y 34 Y 41 Y 42 Y 43 Y 44
and substitute it in (6). Then by computation, we have
Y = Z 1 Z 2 Z 3 Z 4 Z 2 Z 1 Z 4 Z 3 Z 3 Z 4 Z 1 Z 2 Z 4 Z 3 Z 2 Z 1 ,
where
Z 1 = 1 4 Y 11 + Y 22 + Y 33 + Y 44 , Z 2 = 1 4 Y 21 Y 12 + Y 43 Y 34 , Z 3 = 1 4 Y 31 Y 42 Y 13 + Y 24 , Z 4 = 1 4 Y 41 + Y 32 Y 23 Y 14 .
Now, we construct a new quaternion ( R , S )-symmetric matrix X using Z 1 , Z 2 , Z 3 , and Z 4 :
X = Z 1 + Z 2 i + Z 3 j + Z 4 k = 1 4 I m i I m j I m k I m T Y I n i I n j I n k I n   = 1 16 I n i I m j I m k I m T Y + Q m Y Q n T + G m Y G n T + T m Y T n T I n i I n j I n k I n .
It is easy to verify that X τ = Y . By (b) of Lemma 1, we obtain
( A X B ) τ = A τ X τ B τ = A τ Y B τ = C τ .
Thus, X satisfies (1). Moreover, Y is ( R τ , S τ )-symmetric, and we have
( R X S ) τ = R τ X τ S τ = R τ Y S τ = Y = X τ .
It implies that R X S = X , and so X is an ( R , S )-symmetric solution to (1). Therefore, the consistency of (3) implies the consistency of (1). Moreover, any solution Y of (3) can generate an ( R , S )-symmetric solution X of (1).
Next, we show that if (1) is consistent, then (3) is also consistent. Assume (1) has an ( R , S )-symmetric solution X 0 , i.e., R X 0 S = X 0 . By (b) of Lemma 1, we have
A τ X 0 τ B τ = C τ
and
R τ X 0 τ S τ = X 0 τ .
Then we can see that Y 0 = X 0 τ is the ( R τ , S τ )-symmetric solution of (3). Therefore, the original matrix Equation (1) is consistent if and only if its corresponding real matrix Equation (3) is consistent. According to Theorem 2.3 in [7], the real matrix Equation (3) has an ( R τ , S τ )-symmetric solution if and only if the rank equalities in (c) of Theorem 1 hold. Additionally, when (3) is solvable, Y in the form of (5) is the general solution of the real matrix Equation (3), so we can generate our solution X by Equation (4). □
Similar to the proof of Theorem 1, we can obtain the following result for a skew symmetric case.
Theorem 2.
Let A H p × m , B H n × q , C H p × q . Then there are three equivalent statements:
(a) 
The matrix Equation (1) has an ( R , S )-skew symmetric solution X H m × n ;
(b) 
The matrix equation
A τ Y B τ = C τ
has an ( R τ , S τ )-skew symmetric solution Y R 4 m × 4 n ;
(c) 
The following rank equalities hold:
r A 2 A 1 C τ = r A 2 A 1 , r B 1 B 2 C τ = r B 1 B 2 ,
r C τ A 2 B 2 0 = r 0 A 2 B 2 0 , r C τ A 1 B 1 0 = r 0 A 1 B 1 0 .
Furthermore, if the matrix Equation (1) is solvable, then
X = 1 16 I m i I m j I m k I m T Y + Q m Y Q n T + G m Y G n T + T m Y T n T I n i I n j I n k I n
is an ( R , S )-skew symmetric solution to (1), where
Y = P Q 0 Y P V Y Q U 0 U ^ V ^
with
Y P V = Y ˜ P V + S 2 F J U R H T 2 + F A 1 W 1 + W 2 E B 2 ,
Y Q U = Y ˜ Q U + S 1 F J U R H T 1 + F A 2 V 1 + V 2 E B 1 ,
S 1 = I 4 m r 0 , S 2 = 0 I r , T 1 = I k 0 ,
T 2 = 0 I 4 n k , J = A 1 A 2 , H = B 1 B 2 ;
U , V 1 , V 2 , W 1 , W 2 are arbitrary, and Y ˜ P V , Y ˜ Q U are particular solutions to the matrix equation  A 2 Y Q U B 1 + A 1 Y P V B 2 = C τ .

3. Least-Squares ( R , S )-(Skew) Symmetric Solutions

In this section, we derive the least-squares ( R , S )-(skew) symmetric solutions to the quaternion matrix Equation (1). Let M = M 1 + M 2 i + M 3 j + M 4 k H p × m , where M i R p × m . We define
M c τ = M 1 M 2 M 3 M 4 .
We can check that for any N H m × q ,
( M N ) c τ = M τ N c τ .
To simplify the least-squares ( R , S )-(skew) symmetric solution problems, we are going to use the Frobenius norm of M c τ . By direct calculation, we first derive the relation between M c τ F and M | | F :
M | | F = 1 2 M τ | | F = M c τ F .
Let A H p × m , B H n × q , C H p × q . We will find an ( R , S )-(skew) symmetric X = ± R X S H m × n such that
A X B C F = min X 0 H m × n , X 0 = ± R X 0 S A X 0 B C F .
The lemma about the least-squares solutions is given as follows:
Lemma 2.
[22The solutions of the least-squares problem of the complex matrix equation A X = B are
X = A B + ( I A A ) Z ,
in which Z is arbitrary.
Now we recall two important decompositions for an ( R , S )-symmetric (reps. ( R , S )-skew symmetric) matrix over quaternions.
Lemma 3.
[11Let R H m × m , S H n × n be nontrivial involutory matrices with the decomposition
R = E 1 I r 1 0 0 I m r 1 E , S = F 1 I r 2 0 0 I n r 2 F .
Then
(a) 
X H m × n is ( R , S )-symmetric if and only if X can be expressed as
X = E 1 X 1 0 0 X 2 F
where X 1 H r 1 × r 2 , X 2 H ( m r 1 ) × ( n r 2 ) .
(b) 
Y H m × n is ( R , S )-skew symmetric if and only if Y can be expressed as
Y = E 1 0 Y 1 Y 2 0 F
where Y 1 H r 1 × ( n r 2 ) , Y 2 H ( m r 1 ) × r 2 .
Next, we provide the main result of this section.
Theorem 3.
Let A H p × m , B H n × q , C H p × q . Then each least-squares ( R , S )-symmetric solution to the matrix equation A X B = C should be in the form of
X = E 1 X 1 0 0 X 2 F ,
where E and F are given in Lemma 3,
V e c ( X 1 ) = [   I 1 0 1 i I 1 0 1 j I 1 0 1 k I 1 0 1   ] y ,
V e c ( X 2 ) = [   0 2   I 2   0 2   i I 2   0 2   j I 2   0 2   k I 2   ] y ,
y = Pb + (I − PP)z, P = [B1T ⊗ A1, B2T ⊗ A2]τ, b = (Vec(C))cτ,
01 is the zero matrix with the size of r1r2(mr1)(nr2),
02 is the zero matrix with the size of (mr1)(nr2) r1r2,
I1 is the identity matrix with the size of r1r2 × r1r2,
I2 is the identity matrix with the size of k × k,
where k = (mr1)(nr2) × (mr1)(nr2),
z C 4 r 1 r 2 + 4 ( m r 1 ) ( n r 2 )  is arbitrary.
Proof. 
Using the decomposition of an ( R , S )-symmetric matrix X in (a) in Lemma 3, we partition the matrices:
A E 1 = [ A 1 A 2 ] , F B = B 1 B 2 ,
where A 1 H p × r 1 , A 2 H p × ( m r 1 ) , B 1 H r 2 × q , B 2 H ( n r 2 ) × q . Then A X B = C has a least-squares ( R , S )-symmetric solution X if and only if A E 1 X 1 0 0 X 2 F B = C ; that is,
A 1 X 1 B 1 + A 2 X 2 B 2 = C
has a least-squares solution { X 1 , X 2 } . Using the vec operation, we have
( B 1 T A 1 ) V e c ( X 1 ) + ( B 2 T A 2 ) V e c ( X 2 ) = V e c ( C ) ,
which can be rewritten as
[ B 1 T A 1 , B 2 T A 2 ] V e c ( X 1 ) V e c ( X 2 ) = V e c ( C ) .
Taking the real representations on both sides of (9),
[ B 1 T A 1 , B 2 T A 2 ] τ V e c ( X 1 ) V e c ( X 2 ) τ = V e c ( C ) τ .
By Equations (7) and (8),
[ B 1 T A 1 , B 2 T A 2 ] τ V e c ( X 1 ) V e c ( X 2 ) τ ( V e c ( C ) ) τ F = 2 [ B 1 T A 1 , B 2 T A 2 ] τ V e c ( X 1 ) V e c ( X 2 ) c τ ( V e c ( C ) ) c τ F .
Now, we derive the least-squares solution to the linear system
P y = b ,
where
P = [ B 1 T A 1 , B 2 T A 2 ] τ , y = V e c ( X 1 ) V e c ( X 2 ) c τ , b = ( V e c ( C ) ) c τ ,
and Lemma 2 tells us that y = P b + ( I P P ) z is the least-squares solution. If we denote V e c ( X 1 ) = V e c ( X 11 ) + V e c ( X 21 ) i + V e c ( X 31 ) j + V e c ( X 41 ) k and V e c ( X 2 ) = V e c ( X 12 ) + V e c ( X 22 ) i + V e c ( X 32 ) j + V e c ( X 42 ) k , then
V e c ( X 1 ) V e c ( X 2 ) = V e c ( X 11 ) V e c ( X 12 ) + V e c ( X 21 ) V e c ( X 22 ) i + V e c ( X 31 ) V e c ( X 32 ) j + V e c ( X 41 ) V e c ( X 42 ) k .
By the definition of V e c ( X 1 ) V e c ( X 2 ) c τ , we have
y = V e c ( X 1 ) V e c ( X 2 ) c τ = V e c ( X 11 ) V e c ( X 12 ) V e c ( X 21 ) V e c ( X 22 ) V e c ( X 31 ) V e c ( X 32 ) V e c ( X 41 ) V e c ( X 42 ) .
Clearly,
V e c ( X 1 ) = [ I 1 0 1 i I 1 0 1 j I 1 0 1 k I 1 0 1 ] y ,
V e c ( X 2 ) = [ 0 2 I 2 0 2 i I 2 0 2 j I 2 0 2 k I 2 ] y ,
0 1   is the zero matrix with the size of r 1 r 2 × ( m r 1 ) ( n r 2 ) , and 0 2   is the zero matrix with the size of ( m r 1 ) ( n r 2 ) × r 1 r 2 . Then X = E 1 X 1 0 0 X 2 F is our required solutions. □
By using the decomposition in (b) of Lemma 3, we can obtain the general expression of the least-squares ( R , S )-skew symmetric solution as follows.
Theorem 4.
Let A H p × m , B H n × q , C H p × q . Then each least-squares ( R , S )-skew symmetric solution to the matrix equation A X B = C should be in the form of
X = E 1 0 Y 1 Y 2 0 F ,
where E and F are given in Lemma 3, and { A 1 , A 2 } and { B 1 , B 2 } are block matrices determined by A E 1 and F B , respectively.
V e c ( Y 1 ) = [   I 3 0 3 i I 3 0 3 j I 3 0 3 k I 3 0 3   ] y ,
V e c ( Y 2 ) = [   0 4   I 4   0 4   i I 4   0 4   j I 4   0 4   k I 4   ] y ,
y = Pb + (I − PP)z, P = [B2T ⊗ A1, B1T ⊗ A2]τ, b = (Vec(C))cτ,
03 is the zero matrix with the size of r1(nr2) × (mr1)r2,
04 is the zero matrix with the size of (mr1)r2 × r1(nr2),
I3 is the identity matrix with the size of r1(nr2) × r1(nr2),
I4 is the identity matrix with the size of (mr1)r2 × (mr1)r2,
z C 4 n r 1 + m r 2 2 r 1 r 2  is arbitrary.

4. Numerical Example

In this section, a numerical example is provided to illustrate our results. Here, we only consider the ( R , S )-symmetric solution case.
Example 1.
For the given quaternion matrix equation,
A X B = C ,
where
A = 1 + 2 i + 3 j 3 i + j + 2 k 1 + i + 3 k 2 + 3 j + k j + 2 k 2 i + j + 3 k ,
B = 2 i j + k 2 + j + k 3 + i 5 j 4 + 2 i 3 k ,
C = 4 + 10.24 i 0.16 j + 3.4 k 6.92 0.2 i + 9.64 j + 7.44 k 9.52 + 20.96 i + 26.36 j + 3.8 k 0.32 3.68 i + 20.6 j + 17.4 k .
Find the ( R , S ) -symmetric solutions if the matrix equation is solvable, where
R = 1 0 0 0 0.28 0.96 i 0 0.96 i 0.28 , S = 0.28 0.96 k 0.96 k 0.28
are two involution matrices.
Step 1. By Theorem 1, consider the ( R τ , S τ ) -symmetric solution to the corresponding real matrix equation
A τ Y B τ = C τ ,
where A τ R 8 × 12 , B τ , C τ R 8 × 8 , and Y is an unknown matrix.
Step 2. Find P , Q , U , V , such that R τ P = P , R τ Q = Q and S τ U = U , S τ V = V . According to Trench [8], P , Q can be obtained from an orthonormal basis for the eigenspace of R τ associated with λ = 1 , 1 , respectively. Similarly, U , V can be obtained from an orthonormal basis for the eigenspace of S τ associated with λ = 1 , 1 , respectively.
P = 0 0 0 0 1 0 0 0 0.8 0 0 0 0 0 0 0 0 0.6 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0.8 0 0 0 0 0 0 0.6 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0.8 0 0 0 0 0.6 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0.8 0 0 0 0 0 0 0 0 0 0 0.6 0 , Q = 0 0 0 0 0 0.6 0 0 0.8 0 0 0 0 0 0 0 0.6 0 0 0 0 0.8 0 0 0 0 0 0 0 0 0 0.6 0 0 0.8 0 0 0 0 0 0 0 0.6 0 0 0 0 0.8 ,
U = 0 0 0.8 0 0 0 0 0.6 0.8 0 0 0 0 0.6 0 0 0 0.8 0 0 0.6 0 0 0 0 0 0 0.8 0 0 0.6 0 , V = 0 0 0 0.6 0.8 0 0 0 0 0 0.6 0 0 0.8 0 0 0 0.6 0 0 0 0 0.8 0 0.6 0 0 0 0 0 0 0.8 ,
and U ^ , V ^ are calculated from (2).
Step 3. Perform the following decomposition and verify whether (10) has a solution or not.
A τ P Q = A 1 A 2 , A 1 R 8 × 8 , A 2 R 8 × 4 ,
U ^ V ^ B τ = ( B 1 ) T ( B 2 ) T T , B 1 R 4 × 8 , B 2 R 4 × 8 .
By MATLAB, we have
r A 1 A 2 C τ = r A 1 A 2 = 8 , r B 1 B 2 C τ = r B 1 B 2 = 8 ,
r C τ A 1 B 2 0 = r 0 A 1 B 2 0 = 12 , r C τ A 2 B 1 0 = r 0 A 2 B 1 0 = 8 .
Thus, it follows from Theorem 1 that the matrix Equation (10) is consistent.
Step 4. Calculate the ( R τ , S τ )-symmetric solution Y of the matrix Equation (10) by (5):
Y = 0.80 0 0 1.20 1.60 0 0 0.60 1.64 0.48 1.28 0 1.00 0.96 0.36 0.48 0.96 1.00 0.48 0.64 0.48 1.36 0 0.72 0 1.20 0.80 0 0 0.60 1.60 0 1.28 0 1.64 0.48 0.36 0.48 1.00 0.96 0.48 0.64 0.96 1.00 0 0.72 0.48 1.36 1.60 0 0 0.60 0.80 0 0 1.20 1.00 0.96 0.36 0.48 1.64 0.48 1.28 0 0.48 1.36 0 0.72 0.96 1.00 0.48 0.64 0 0.60 1.60 0 0 1.20 0.80 0 0.36 0.48 1.00 0.96 1.28 0 1.64 0.48 0 0.72 0.48 1.36 0.48 0.64 0.96 1.00 .
Step 5. Generate X from Y by the formula (4):
X = 0.80 0 1.64 0.48 0.96 1.00 + 0 1.20 1.28 0 0.48 0.64 i + 1.60 0 1.00 0.96 0.48 1.36 j + 0 0.60 0.36 0.48 0 0.72 k .
Step 6. Verify whether X is an ( R , S ) -symmetric solution to A X B = C or not. By direct computation, X satisfies
R X S = X , A X B = C .
Thus, X is our required solution.

Author Contributions

Conceptualization, R.L.; methodology, R.L.; validation, S.L.; formal data curation, X.L.; writing—original draft preparation, R.L.; writing—review and editing, R.L., X.L., S.L. and Y.Z.; supervision, Y.Z.; funding acquisition, X.L., S.L. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the Macao Science and Technology Development Fund (No. 0013/2021/ITP); the grants from the National Natural Science Foundation of China (12371023, 12271338, 12001259) and the Natural Sciences and Engineering Research Council of Canada (NSERC) (RGPIN 2020-06746); the Joint Research and Development Fund ofWuyi University, Hong Kong and Macao (2019WGALH20); the MUST Faculty Research Grants (FRG-22-073-FIE); the Science Foundation of Fujian Province (2020J01846); and the Research Foundation of Minjiang University for the Introduction of Talents (MJY17006).

Data Availability Statement

The data supporting the findings of this study cannot be made publicly available due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chen, Y.; Jia, Z.G.; Peng, Y.; Zhang, D.A. New structure-preserving quaternion QR decomposition method for color image blind watermarking. Signal Process. 2021, 185, 108088. [Google Scholar] [CrossRef]
  2. Gaudet, C.J.; Maida, A.S. Deep quaternion networks. In Proceedings of the 2018 International Joint Conference on Neural Networks, Rio de Janeiro, Brazil, 8–13 July 2018; pp. 1–8. [Google Scholar]
  3. He, Z.H.; Wang, Q.W. A real quaternion matrix equation with applications. Linear Multilinear Algebra 2013, 61, 725–740. [Google Scholar] [CrossRef]
  4. Mandic, D.P.; Jahanchahi, C.; Took, C.C. A quaternion gradient operator and its applications. IEEE Signal Process. Lett. 2011, 18, 47–50. [Google Scholar] [CrossRef]
  5. Yang, Y. Spacecraft attitude determination and control: Quaternion based method. Annu. Rev. Control 2012, 36, 198–219. [Google Scholar]
  6. Zhu, X.; Xu, Y.; Xu, H.; Chen, C. Quaternion convolutional neural networks. Proc. Eur. Conf. Comput. Vis. 2018, 7, 631–647. [Google Scholar]
  7. Dehghan, M.; Hajarian, M. Matrix equations over (R, S)-symmetric and (R, S)-skew symmetric matrices. Comput. Math. Appl. 2010, 59, 3583–3594. [Google Scholar] [CrossRef]
  8. Trench, W.F. Minimization problems for (R, S)-symmetric and (R, S)-skew symmetric matrices. Linear Algebra Appl. 2014, 389, 23–31. [Google Scholar]
  9. Chen, H.C. Generalized Reflexive Matrices: Special Properties and Applications. SIAM J. Matrix Anal. Appl. 1998, 19, 140–153. [Google Scholar] [CrossRef]
  10. Lv, C.Q.; Ma, C.F. The iterative algorithm for solving a class of generalized coupled Sylvester-transpose equations over centrosymmetric or anti-centrosymmetric matrix. Int. J. Comput. Math. 2019, 96, 1576–1594. [Google Scholar]
  11. Wang, Q.W.; Van der Woude, J.W.; Chang, H.X. A system of real quaternion matrix equations with applications. Linear Algebra Appl. 2009, 431, 2291–2303. [Google Scholar] [CrossRef]
  12. Xu, W.R.; Chen, G.L. Inverse problems for (R, S)-symmetric matrices in structural dynamic model updating. Comput. Math. Appl. 2016, 71, 1074–1088. [Google Scholar] [CrossRef]
  13. Xie, Y.; Huang, N.; Ma, C. Iterative method to solve the generalized coupled Sylvester-transpose linear matrix equations over reflexive or anti-reflexive matrix. Comput. Math. Appl. 2014, 67, 2071–2084. [Google Scholar]
  14. Yuan, Y.; Dai, H. Generalized reflexive solutions of the matrix equation AXB = D and an associated optimal approximation problem. Comput. Math. Appl. 2008, 56, 1643–1649. [Google Scholar]
  15. Zhang, Q.; Wang, Q.W. The (P, Q)-(skew) symmetric extremal rank solutions to a system of quaternion matrix equations. Appl. Math. Comput. 2011, 217, 9286–9296. [Google Scholar] [CrossRef]
  16. Peng, Z.Y.; Hu, X.Y. The reflexive and anti-reflexive solutions of the matrix equation AX = B. Linear Algebra Appl. 2003, 375, 147–155. [Google Scholar] [CrossRef]
  17. Wang, Q.W.; Chang, H.X.; Lin, C.Y. P-(skew) symmetric common solutions to a pair of quaternion matrix equations. Appl. Math. Comput. 2008, 195, 721–732. [Google Scholar] [CrossRef]
  18. Zhou, F.Z.; Zhang, L.; Hu, X.Y. Least-square solutions for inverse problems of centrosymmetric matrices. Comput. Math. Appl. 2003, 45, 1581–1589. [Google Scholar] [CrossRef]
  19. Jiang, T.S. Algebraic methods for diagonalization of a quaternion matrix in quaternionic quantum theory. J. Math. Phys. 2005, 46, 783–786. [Google Scholar] [CrossRef]
  20. Liu, X.; Wang, Q.W.; Zhang, Y. Consistency of quaternion matrix equations AXXB = C,XAXB=C. Electron. J. Linear Algebra 2019, 35, 394–407. [Google Scholar] [CrossRef]
  21. Rodman, L. Topics in Quaternion Linear Algebra. In Princeton Series in Applied Mathematics; Princeton University Press: Princeton, NJ, USA, 2014. [Google Scholar]
  22. Ben-Israel, A.; Greville, T.N.E. Generalized Inverses: Theory and Applications; Wiley: New York, NY, USA, 1974. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liao, R.; Liu, X.; Long, S.; Zhang, Y. (R, S)-(Skew) Symmetric Solutions to Matrix Equation AXB = C over Quaternions. Mathematics 2024, 12, 323. https://doi.org/10.3390/math12020323

AMA Style

Liao R, Liu X, Long S, Zhang Y. (R, S)-(Skew) Symmetric Solutions to Matrix Equation AXB = C over Quaternions. Mathematics. 2024; 12(2):323. https://doi.org/10.3390/math12020323

Chicago/Turabian Style

Liao, Ruopeng, Xin Liu, Sujuan Long, and Yang Zhang. 2024. "(R, S)-(Skew) Symmetric Solutions to Matrix Equation AXB = C over Quaternions" Mathematics 12, no. 2: 323. https://doi.org/10.3390/math12020323

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop